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We consider random transverse-field Ising spin chains and study the magnetization and the energy- 
density profiles by numerically exact calculations in rather large finite systems [L < 128). Using 
different boundary conditions (free, fixed and mixed) the numerical data collapse to scaling functions, 
which are very accurately described by simple analytic expressions. The average magnetization 
profiles satisfy the Fisher-de Gennes scaling conjecture and the corresponding scaling functions are 
indistinguishable from those predicted by conformal invariance. 
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Every experimental system is geometrically con- 
strained and therefore has a surface, for which reason 
we have to discriminate between so called bulk and sur- 
face properties. This is justified as long as the correlation 
length is much smaller than the system size. However, at 
a critical point it is more appropriate to describe the posi- 
tion dependent physical properties of the system by den- 
sity profiles rather than bulk and/or surface observables. 
For a number of universality classes much is known about 
this spatially inhomogeneous behaviorH, in particular in 
two dimensions, where conformal invariaxice provides a 
powerful tool to study various geometriesB. 

Not much is known about this issue for quantum sys- 
tems with quenched (i.e. time independent) disorder. 
Here one is confronted with a possible quantum phase 
transition, i.e. a zero temperature transition that is trig- 
gered by quantum rather than thermal fluctuations, as 
for instance in random transverse Ising modelsu lI. Their 
bulk properties have been studied quite extensively by 
now. The aim of the present letter is to investigate for 
the first time the above mentioned density profiles in a 
geometrically constrained disordered system at a quan- 
tum phase transition. In particular we study numerically 
the random transverse field Ising chain and propose ana- 
lytic expressions of the magnetization and energy density 
profile for various boundary conditions (b.c). 

In a critical system confined between two parallel 
plates, which are a large but finite distance L apart, the 
local densities {^{r)) such as the order parameter (mag- 
netization) or the energy density vary with the distance 
I from one of the plates as a smooth function oil/L. 
cording to the scaling theory by Fisher and de Genned 
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where is the scaling dimension of the operator <&, while 
ab denotes the boundary conditions at the two plates. 
The scaling function in (|l|) has the asymptotic behavior: 
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where the exponent in the fiE£lj-|Correction term was con- 
firmed by different jnethodaJcl. It_has been shown by 
Burkhardt and XueE3 and by CardyQ that the Bab coeffi- 
cients in (^ and the Aab finite size correction coefficients 
of the free energy as AabL~'^'^^ are related to each other: 
their ratio is universal and independent of the form of 
the b.c. 

Having the same type of b.c. at both plates the pro- 
file {^{l))aa = L^^^ faa{ll L) is rcficction symmetric 
faa{v) = /aa(l " v) and accordiug to eqs.(|l]) and (|) 

lim„^o /aa(i') . Conseauently [/aa(w)]~^^''* can 

be expanded in a Fourier seriesllll and the profile is given 
by: 
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The Fourier expansion in (^ has different convergence 
properties in two- and three-dimensions due to the differ- 
ent parity of the correction term in (^. While in three- 
dimensions infinite terms are needed to reproduce the 
Fisher- de Gennes scaling result in (|^) , in two-dimensions 
one expect to obtain satisfactory accuracy by the first few 
terms of the expansion. Indeed for conformally invariant 
two-dimensional models only the first term inJjhe Fourier 
series in (^) gives non-vanishing contributioncl: 
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Conformal invariance can be used further to predict the 
density profiles with geiijsral b.c. In two-dimensions the 
profiles are in the forn£3 
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where the scaling function Gabil/L) depends on the uni- 
versality class of the model and on the type of the b.c. 
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For the Ising model the magnetization profile with free- 
fixed boundary condition the scaling function is predicted 



Gf+ = B 
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where = 1/2 is the scaling dimension of the surface 
magnetization operator. Similar result is obtained for the 
Q < 4 state Potts modelllj with the appropriate surface 
scaling dimension in (j^). 

In the present Letter we consider the random trans- 
verse filed Ising chain 



H = - 
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Here the J/ exchange couplings and the hi transverse- 
fields are independent random variables with distribu- 
tions 7r(J) and p{li), respectively and the erf, erf are 
Pauli matrices at_site I. This Hamiltonian is the extreme 
anisotropic limitll3 of the layered two-diapmsional Ising 
model as introduced by McCoy and Wut3'E3. 

The critical behavior of the random transverse- 
field Ising rSpiii-, chain in (|^) haa_.been investigated 
analyticallytj'EjQ and numericallyll30 in several papers. 
Depending on the strength of the average value of the 
transverse-field the system has two phases, which are 
separated,by a second order phase transition point lo- 
cated atlla 5 = In J — In ft, = 0. Due to a broad dis- 
tribution of various physical quantities the typical and 
average quantities of the system are generally different. 
The scaling dimensions of,the averaged rrvagnetization are 
Xm = (3-75)74 « 0.191Q and a;^ = 1/2^. The model is 
anisotropic at the critical point, the dynamical exponent 
is z = oo. More precisely the characteristic length scale 
^ and the corresponding time scale t are related through: 
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thus the model is not conformally invariant and predic- 
tions in eqs. (H-^) are not expected to be valid. 

In the following we briefly describe how the den- 
sity profiles were calculated. The local magnetiza- 
tion mi is obtained from the asymptotic behavior of 
the (imaginary) time-time correlation function Gi(t) = 
{af{j)afm = m^m? ^M-riE, - E^)] where 
(0| and (i| denote the ground state and the z-th excited 
state with energies Eq and Ei^ respectively. In the low 
temperature (strong coupling) phase Ei is asymptotically 
degenerate with the ground state, thus the sum is domi- 
nated by the first term. In the large r limit &(''") = 
therefore the local magnetization is given byllS 



I ' 



mi 



(ikflo) 



(9) 



The energy-density profile is given by the ground state 
expectation value e; = (0|crf |0). Since e; contains a non- 
singular contribution the scaling behavior of the energy- 
density is more convenient to deduce from the asymptotic 



form of the connected time-time correlation function of 
the energy-density operator erf. Similarly to the order- 
parameter the singular energy density e; is given by a 
matrix-element: 
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where (e| denotes the first excited state, which has non- 
vanishing matrix-element with the ground state. 

To calculate the matrix-elements in eqs.|-(P) and (0) 
we first, following Lieb et a£j and PfeutyE37 transform 
H into a free-fermion model. For the fixed and free b.c. 
we study in this letter we found it npst convenient to 
choose the representation described inO, which necessi- 
tates only the diagonalization of an 2L x 2L-tridiagonal 
matrix. From the corresponding eigenvectors one obtains 
the local magnetization (^ and local energy density (|l0| ) 
as described inE2l. 

The critical properties of random Ising chains are ex- 
pected to be independent of the details of the distribu- 
tions of the couplings and/or fields. In this Letter we 
consider two different cases: the binary distribution 

7r(J) = i<5(J-A) + i,5(J-A-i); h = ho , (11) 

i.e. p{h) — 5(h — ho), and the uniform distribution: 

7r(J) = 61(1 - J)e{J) ■ p{h) = 0{ho - h)9{h) (12) 

In both cases the critical point is at /ig = 1- AH numerical 
data which we present below are averaged over 50,000 
samples and the resulting statistical error is much smaller 
than the size of the symbols used in the plots. Disorder- 
averaged quantities are denoted by the brackets [. . .]av 

First we study the magnetization profile of the system 
with fixed b.c. at both ends of the chain. The finite size 
results on the pure model, which are shown in the inset 
of Fig. 1 are in complete agreement with the conformal 
prediction in (^. The profile for the random chain is 
shown in Fig. 1. From the scaling plot one can see that 
the Fisher-de Gennes scaling result in (|l|) is well satisfied 
with the conjectured value of the decay exponent Xm — 
(ijv — .191. Note that we do not use Xm (as well as later 
x^) as fit paramters but fix them to the theoretically 
predicted values cited above. The only fit parameter is 
the non-universal prefactor A'yh (^) . Obviously, one can 
very accurately describe the finite-size data in the whole 
profile with the first term of the Fourier expansion in (|^). 
The corrections to the conformal result in (|^) are indeed 
negligible. 

Next we turn to study the magnetization profiles with 
free-fixed b.c. As seen on the inset of Fig. 2 the finite- 
lattice results on the pure model perfectly coincide with 
the conformal prediction in (^ . Results for random mod- 
els are shown in Fig. 2. As one can see the numerical data 
collapse to a scaling function, which can be very accu- 
rately described by a function of the form in (^) with the 
exponents: Xm — -191 and ~ (i'^ /v = 1/2 (again the 
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FIG. 1. Scaling plot of the magnetization profile mL(Z) (y) 
with fixed b.c. at both ends. We have shifted the site index by 
half a lattice constant and denote I' = 1 — 0.5. The pure case 
is depicted in the inset, for which the scaling function is given 
by (|) with x^'''' = 1/8 and G++ = const. The main figure 
shows the result for the binary distribution (|l^) with A = 4. 
Other values of A as well as the uniform distribution yield the 
same quality for the data collapse, with different values for 
the non-universal prefactors but identical scaling function (^) 
with s™"^"™ = P/v ^ 0.191 and G++ = const. 



only fit-parameter is the non-universal prefactor A ■ B 
from According to Fig. 2 the corrections to the 

conformal result seem here also to be negligible. 




FIG. 2. The same as in fig. 1 with fixed b.c. on the right 
end of the chain and free b.c. on the left end. The magnetiza- 
tion profile is given by ^ and with a;^ = 1/2 for the pure 
and the random case. The data shown in the main figure are 
for the uniform distribution ([l2|). 

Due to symmetry the magnetization in a finite system 
with free (non-symmetry breaking) b.c. is zero. However, 
introducing an infinitesimal symmetry breaking field ii_ 
one can obtain the magnetization by Yang's methodcj 
with the identical result as in (^). Then, as h ^ the 
matrix-element in can be considered to define the 
local magnetization in a finite system in a time scale 



T <C where t is the relaxation time. The profile of 
the matrix-element in (|^) can be predicted by conformal 
invarianceo. For a general local operator ^{l) the scaling 
form in the strip geometry is given bynil: 



(13) 



where 



denotes the surface scaling dimension of 



$. This expression satisfies the known scaling limits 
(0|$(1)|$) - L-^i and (0|$(L/2)|$) - L^^"* at the sur- 
face and in the bulk, respectively. For non-conformally 
invariant systems ( [l^ ) represents the first leading term 
of a Fourier-expansion, as in (p|) and (0). 
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FIG. 3. The same as in fig. 1-2 with free b.c. on both ends 
of the chain. The data for the random case are for a binary 
distribution with A = 2. The magnetization profile is given by 
(p^. The data shown in the main figure are for the uniform 
distribution (p^. 

Numerical results on the magnetization profiles with 
free boundary conditions are shown on Fig. 3. Again the 
finite-size results on the pure Ising model are in complete 
coincidence with the conformal prediction in (|l3|). For 
the random case the numerical data collapse to a scaling 
curve, which is very accurately described by the confor- 
mal expression in ( |l3| ) with the exponents Xm — -191 and 
= 1/2- Thus again the non-conformal corrections are 
very small. 

Finally, we discuss the singular part of the energy den- 
sity profile and study the energy density matrix-element 
in (p^. For the pure model one can easily evaluate e(Z), 
which yields in the scaling limit (/ ^ 1, i ^ 1): 



eil) = — sinTT— 
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with 



This corresponds to the conformal result in eg die 
Xe = 1 and = 2. 

In a quantum system the bulk energy-density is pro- 
portional to the inverse relaxation time: e ~ t^^. In 
the random transverse Ising chain the scaling is anoma- 
lous as indicated in (||), therefore the appropriate scaHng 
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combination is L~^/^ Ine(Z) instead of Le{iy^^ if z would 
be finite. In the following we study the typical energy 
density [lne(Z)]^y, which after multiplication with L^^/"^ 
yields a universal scaling function. The finite size data 
for the random case very well satisfy the relation: 



[lne(/)Ui:i/2 = A„ + Ai(^sin7r-^' 



(15) 



We note that this expression can also be considered as the 
leading part of a Fourier expansion, where the correction 
terms are again very small. 




128 

-3.03+1.33*sin(7ix)"-^ 



0.2 



0.4 



0.6 



0.8 



I'/L 



FIG. 4. The energy profile [l) ( [lo| ) with free b.c. on both 
sides. The pure case is depicted in the inset, for which the 
scaling function is given by (|lj). The main figure shows the 
result for the binary distribution (^) with A = 4. Here the 
scaling function is well described by (|l5|). We note that the 
approach to the asymptotic scaling limit seems to be much 
slower than for the magnetization profiles. 

To summarize we have investigated the density profiles 
of random trans verse- field Ising spin chains. The numer- 
ical data on rather large systems L < 128 follow scaling 
plots and the scaling functions can be described very ac- 
curately by analytical expressions, which are derived for 
conformally invariant systems. Since our system is not 
conformally invariant there are presumably corrections. 
These are, however, very small, certainly smaller than 
the error in our present numerical calculation. 

Generally the non-conformal corrections to the den- 
sity profiles are not small. As an example we mention 
the two-dimensional aperiodically layered Ising modeled, 
which is somewhat related to our problem. When the 
aperiodically modulated couplings of the model repre- 
sent a marginal perturbation the system is desciihed by 
a coupling dependent dynamical exponent z > lEll, thus 
the system is not conformally invariant. Although the 
aperiodic model looks similarly to our random problem 
its density profiles are completely different from the con- 
formal resultsE3. One could speculate about the existence 
of some hidden symmetry which explains the coincidence 
of the density profiles of the random transverse-field Ising 
chain with the conformal result. 
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